Hysteresis and hierarchies: dynamics of 
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ABSTRACT: We use the zero-temperature random-field Ising model to study hys- 
O' 

2 I teretic behavior at first-order phase transitions. Sweeping the external field through zero, 

Q\ , the model exhibits hysteresis, the return-point memory effect, and avalanche fiuctuations. 
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There is a critical value of disorder at which a jump in the magnetization (corresponding 



F^ ■ to an infinite avalanche) first occurs. We study the universal behavior at this critical point 

a: 

O ■ using mean-field theory, and also present preliminary results of numerical simulations in 



three dimensions. 
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First-order phase transitions have always been the weak sibhng of critical phenomena 
in statistical mechanics. At some critical temperature Tc, the idealized thermodynamic 
equilibrium properties of a homogeneous material abruptly shift from (say) liquid to gas, 
with precursor fluctuations almost entirely absent. The classic model for a first order 
transition is the Ising model in an external field H at T < Tq. as H passes through zero, 
the equilibrium magnetization reverses abruptly. 

There is an amazing contrast, though, with real first-order transitions as studied by 
materials scientists and metallurgists (figure 1). A solid material which transforms from 
one crystalline or magnetic form to another under the infiuence of temperature, external 
stress, or applied field often has no sharp transition at all. Hysteresis becomes the key 
phenomenon. Also, there are the experimental Barkhausen and return-point memory 
effects, which are explained here in terms of avalanches and a hierarchy of metastable 
states. 

The ingredient we will add to the idealized physics picture is disorder. Adding a 
random field fi at each site of the Ising model 

l-L = -^JijS^Sj-^f^s^ + Hs^ (1) 

ij i 

allows us to study the effect of disorder in a first order phase transition. We study this 
random-field Ising model at zero temperature. As the external field H is changed, each 
spin will flip when the direction of its total local field 

Ft = Yl -^^J^J + h + H (2) 

3 

changes. The system transforms from negative to positive magnetization as the field is 

swept upwards. Related approaches ' have been useful for studying the Barkhausen effect, 
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but they have not discussed the return-point memory and did not attempt to vary the 
disorder and field to reach a critical point. 

The dynamics of the approach to thermal equilibrium in this model has been well 
studied . However, here we study the athermal dynamical response to an external field, 
where the transition comprises a series of cluster flips. This model is applicable to many 
experimental systems where the elementary domains have barriers to flipping large enough 
that thermal activation can be ignored. For example, martensitic transitions come in two 
varieties, called isothermal and athermal: the athermal martensites show no transformation 
if the temperature and external strain are held constant, but transform with a crackling 
noise as these parameters are changed. Similarly, most useful magnetic memory devices 
by design do not come to thermal equilibrium! 

Simulating the system described in equation (1) at zero temperature yields behavior 
pleasantly familiar to the experimentalist. Figure 1 shows the hysteresis loop for a three- 
dimensional 30 system with nearest-neighbor bond strength J = 1 and random flelds /j 
with a Gaussian probability distribution with a half-width R. The outer loop shows the 
external field H being swept from a large negative value (saturating all spins to -1) to a 
large positive value and back. 

The inner loops in figure 1 illustrate the return-point memory effect, seen in some but 
by no means all first order transitions . (Some first-order transitions exhibit a drift in 
their hysteretic behavior. ) If the field H is made to backtrack from Hj^ to Hq, when it 
returns to Hj^ the system returns precisely to the same state from which it left the outer 
loop, it remembers the former state. The slope of the M{H) curve has a nonanalyticity as 
it rejoins the outer loop: when the field H{t) passes through the previous local maximum 
Hj^, new spins begin to fiip, leading to the apparent slope discontinuity at point B. This 
same memory effect extends to subcycles within cycles (and so on); the state of the system 
can remember an entire hierarchy of turning points in its past external field H{t). 
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This memory effect is vividly illustrated by experiments measuring the avalanche fluc- 
tuations during the transformation. The nucleation of the individual domains in marten- 
sites can produce observable pulses in acoustic emission and latent heat. Under a cycle 
like that between Hj^ and H(j in figure 1, the acoustic emission is resolved into hundreds 
of individual pulses each lasting some /usec, of varying height, forming a "fingerprint" of 
the cycle. On repeating the cycle, the avalanche structure is precisely reproduced! 

We can explain the return-point memory effect theoretically using Alan Middleton's 
"No Passing" rule, introduced in the study of sliding charge-density waves. Consider the 
natural partial ordering of the states: a state s = {si,...,S7v}>r = {^'i, • • • , rj^} if Si > ri 
for each site i in the system. This ordering is not very discriminating: most arbitrary pairs 
of states will not have any definite relationship. For example, in figure 1, the states a and 
b are unrelated: despite the fact that the net magnetization of a is smaller than that of 
b, there is likely at least one spin that has fiipped on the way from a to B which has 
not flipped back down on the way down to b. This partial ordering becomes important 
because it is preserved by the dynamics: 

No Passing: (Middleton) Let a system s(t) be evolved under the field Hs{t), and similarly 
r(t) evolved under Hrit). Suppose the initial configurations satisfy s(0) > r(0) and the 
fields Hgif) > Hr{t). Then it will remain true that s(t) > r(t) at all times t. 

Proof: Suppose the contrary. Then there must be some first time that a spin j in r is 
going to "pass" the corresponding spin in s, i.e. rj{t) > Sj{t). At that time, the local field 
F^ must be larger than F? But this cannot be, because the neighbors of tj in state r at 
that time are no less negative than those of Sj and the external field Hr < Hg. 

Having established the No Passing property of the system, we make the additional 
assumption that the system dynamics is adiabatic: the field changes slowly enough that 
if we start in some state A any monotonic path from field Hj{ to Hj^ will cause the state 
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to evolve in the same way, and into the same final state B. For this model and for other 
systems with partial ordering, No Passing, and adiabaticity, we can now prove: 

Return— Point Memory: Suppose a system s(0) is evolved under field H{t), where 
H{0) < H{t) < H{T) for < t < T, with H{t) not necessarily monotonic. Then the 
final state of the system depends only on H{T), and is independent of the time T or the 
history H{t). In particular, a system coming back to a previous extremal field will return 
to exactly the same state, provided that the field remains within these bounds. 

Proof: Consider the fields H-^{^{t) = min^/>^ ii/'(i') and -ffmax(i) = iiiax^/<^ii/'(t'). These 
fields irregularly but monotonically rise from H{Q) to H{T). If Sniin evolves under i^niin 
and similarly Smax evolves under -ffniax, then since H^{^{t) < H{t) < Hjjiax{i) the No 
Passing rule implies Sjj^jj^(t) < s(t) < Sniax(^)- But by the adiabaticity assumption, all 
monotonic paths lead to the same final state Sjnin(T) = Smaxl^"), so s(T) is independent 
of path. 

Historically, ours is not the first model to exhibit hysteresis and return-point memory. 
These phenomena are often described by Preisach models , where the system is decom- 
posed into independent elementary hysteresis domains, each with an upper and lower crit- 
ical field for flipping. Preisach models have the three properties needed for return-point 
memory, and indeed can be thought of as zero-dimensional variants of our model. They are 
successful because they demonstrate return-point memory, and because the distribution 
of domains can be varied in order to flt experimental hysteresis loops. The hypothesis of 
independent domains, of course, is an idealization, since there is generally some non-trivial 
coupling between components of the system. As an illustration of the collective behavior 
missing in the Preisach model, we unearth critical fluctuations and universality buried in 
the dynamics of the interacting system. 

Physically, if the nearest-neighbor coupling is substantial compared to the randomness, 
i.e. bond strength J^R, then the system will look like what one expects of a clean first 
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order transition: the first sufficiently large region to nucleate will push most of the rest 
of the sample over the brink. In this small disorder regime, there will be an "infinite 
avalanche" : some spin fiips, triggering its neighbors to fiip (and so on) until a finite fraction 
of the sample is transformed, causing a jump in the magnetization density. This infinite 
avalanche will perhaps be surrounded by precursors and aftershocks. On the other hand, 
in the large disorder regime, where J<^R the spins will essentially fiip on their own, each 
spin hipping as the external field crosses its local random field. Those avalanches which 
occur will be small: there will be no diverging correlation lengths. 

Consider the "phase diagram" of our model as we vary disorder R and field H (at 
T = 0), starting from the state with all spins down (figure 2). There must be a strength 
of disorder R^ above which an infinite avalanche never happens. In contrast, for weaker 
disorder, there will be a field H^{R) at which (on raising the field from the down state) 
an infinite avalanche occurs. Both experimentally and in our initial simulations in three 
dimensions, the transition at H^(R) seems to be abrupt for R < Rq. (We'd call it first- 
order, except the avalanche is supposed to be mediating a first-order transition — the 
language is failing us.) On the other hand, as one approaches the critical field H^{Rc) at 
the critical disorder R^ the transition appears to be continuous: the magnetization M{H) 
has a power-law singularity, and there are avalanches of all sizes. As one approaches this 
"endpoint" at {R^ H'^{Rc)) in the (i?, H) plane, we find diverging correlation lengths and 
presumably universal critical behavior. 

We can solve for these critical properties exactly within mean-field theory. Suppose 
every spin in equation (1) is coupled to all A^ other spins with coupling J/N. The effective 
field acting on a site is JM + f^ + H, where M = ^^ s^/N is the average magnetization. 
Spins with /^ < —JM — H will point down, the rest will point up. Thus the magnetization 

6 



M{H) is given implicitly by the expression 

j'-JM{H)-H 
M{H) = 1 - 2 / p{f) df (3) 

where p{f) is the probability distribution for the random field. This equation has a single- 
valued solution unless R < Re (which in the case of a Gaussian distribution corresponds 
to p{0) > 1/2 J), at which point hysteresis and an infinite avalanche begin. Near the 
endpoint, the jump in the magnetization AM due to the avalanche scales as r^, where 
r = [Rq—R)/ Re and [3 = 1/2. As one varies both r and the reduced field h = [H—W^^Rc)), 
the magnetization scales as M(h^r) ~ \rY M.±{h/\rY )^ where ± refers to the sign of r. 
In mean-field theory 5 = 3 and M.± is given by the smallest real root g±{y) of the cubic 
equation g^ ^ ^g - -^j^y = 0. 

Unfortunately, the mean field theory predicts unphysical behavior in two ways. First, 
there is no hysteresis apart from the infinite avalanche. This is an artifact of the Ising 
mean-field theory, and is not observed in finite dimensions nor in an otherwise equivalent 
soft-spin mean-field model. Second, the approach to the infinite avalanche upon varying 
H for R < Re is continuous in mean-field theory. However, in three dimensions we 
numerically observe a discontinuous transition, although fiuctuation effects do seem to be 
large. 

More interesting is the avalanche size distribution near the critical point, (inset, fig- 
ure 2). We can solve exactly for the probability D{s,t) of having an avalanche of size s, 
where t = 2Jp{—JM — H) — 1 measures the distance to the infinite avalanche at p = 1/2 J. 
To have an avalanche of size s triggered by a spin with random field /, you must have 
precisely s — 1 spins with random fields in the range {/, / + 2Js/N}. The probability for 
this to happen is given by the Poisson distribution. In addition, they must be arranged 
so that the first spin triggers the rest. This occurs with probability precisely 1/s, which 
we can see by putting periodic boundary conditions on the interval {/, / + 2Js/N} and 
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noting that there would be exactly one spin of the s which will trigger the rest as a single 
avalanche. This leads to the avalanche size distribution 

D{s,t) = j^_{t + ir-'e-<'+'\ (4) 

To put this in a scaling form, we must first express t as a function of r and h: 



t r^ r 



1 T ffi'±(^/|^h) . Using some simple expansions and Stirling's formula, we can 



then write D in the scaling form 

D{s,r,h) ^ s-^V±{s/\r\-^/'',h/\rf^) (5) 

where our mean-field calculation gives r = 3/2, a = 1/2, and the universal scaling function 

V±{x, y) = ^e-^[lTi5±(y)'] Vs. (6) 

V 27r 



As usual, we expect the critical exponents /3, r, 5, and a and the scaling functions M. 
and V to be independent of many details of the system (and thus the same for theoretical 
models and real materials), but to depend on dimension, the range of the interaction, 
and on the symmetries of the order parameter. We intend to extract critical exponents 
both numerically and via an e-expansion (the upper critical dimension for the closely 
related charge-density-wave depinning problem is four ), study finite temperatures, and 
introduce frustration representing the dipole fields in magnets and the elastic strain in 
martensites. 

Does the critical behavior predicted here exist in the real world? An FeNi alloy showed 
a crossover from athermal to burst (infinite avalanche) behavior as the grain size was 
varied. Grain boundaries are not random fields, but we expect critical fluctuations and 
scaling where the bursting first occurs. The distribution of avalanches in magnetic systems 
has been studied, and some preliminary fits to power laws have been made. Avalanches 
and hierarchies have been implicit in the literature since the 1920's; power laws and critical 
scaling are the tools needed for understanding the collective behavior being studied now. 
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Figure Captions 

Figure 1: Hysteresis Loop Showing Return— Point Memory. Shown is the mag- 
netization as a function of external field for a 30 system with disorder R = 3.5 J. Note 
that the system returns to the original curve at exactly the same state B that it left, that 
the returning curve has an apparent slope discontinuity at B, and that both effects also 
happen for the internal subloop. Thus a state can have a whole hierarchy of parent states 
(mothers at increasing fields and fathers at decreasing fields), which are seen as kinks in 
the corresponding branch of the H{M) curve. 

Figure 2: Varying the Disorder. Three H{M) curves for different levels of disorder, for 
a 60 system. Our current estimate of the critical disorder is Re = 2.21J (we set J = 1). 
At R = 2 < Re, there is an infinite avalanche which seems to appear abruptly. For R = 
2.6 > Re the dynamics is macroscopically smooth, although of course microscopically it is 
a sequence of sizable avalanches. At R = 2.3, near the critical level of disorder, extremely 
large events become common. Inset: Log-Log Plot of the avalanche-size distribution D{s) 
vs. avalanche size s for the 60 system at R = 2.3 for 1.3 < H < 1.4, averaged over 20 
systems. Here D{s) ~ s~ , compared to the mean-field exponent r of 3/2. We expected 
to see a cutoff in the avalanche sizes: presumably we will see it farther from the endpoint 
or for larger system sizes. 
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